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ABSTRACT 


The inherent brittleness of glass invariably leads to a 
large variability in strength data and a time dependence in 
strength. Loading rate plays a large role in strength values. 
Glass is found to be weaker when supporting loads over long 
periods of time as compared to glass which undergoes rapid 
loading. These properties complicate the structural design 
allowables for the utilization of glass components in an 
application such as the AXAF. 

This report describes the test methodology to obtain 
parameters which can be used to predict the reliability and 
lifetimes of Zerodur glass-ceramic which is to be used for the 
mirrors in the AXAF. 
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INTRODUCTION 


The aim of this study is to provide the necessary data to 
predict the reliability and lifetime of Zerodur Glass-Ceramic 
which will be utilized as the material for the grazing incidence 
mirrors in the Advanced X-ray Astrophysics Facility (AXAF) . 

Zerodur is a alumino-silica glass-ceramic containing 70 to 
78% by weight of high-temperature quartz as the crystalline 
material. The mean crystal size of the crystals is in the region 
of 50 to 55 NM. The crystalline phase has a negative expansion 
coefficient, while the vitreous phase has a positive expansion 
coefficient. This pratically balances out linear thermal 
expansion leading to long term stability even with temperature 
fluctuations. Zerodur, is therefore an ideal candidate for AXAFS 
and other applications such as mirror substrates mounts where 
changes in the material due to temperature variations could 
impair the quality of observations. The composition and 
mechanical, thermal and optical properties have been reported by 
Schott Glass Technologies, producers of Zerodur. However, the 
question of subcritical crack growth in this material has not 
been explored in detail. 

The design of the stresses in the AXAF mirrors is dictated 
by the ultimate strength of the material and a nominal value of 
1450 psi (MOR) is recommended by the manufacturer. This is a 
minimum value or worst case value. In glassy materials, the 
modulus of rupture depends upon primarily the surface condition 
of the material. In other words, a surface with a large number 
of flaws will have a lower strength than one with few on no 
flaws. The reduction in strength is caused by surface flaws 
which concentrate applied stresses leading to failure at much 
lower than theoretical applied loads. These flaws can be 
introduced by machining or normal handling. Thus one method to 
increase strength is to polish the surface thereby removing 
surface flaws. Schott Glass reports an MOR of approximately 
11,000 psi with a failure probability of 0.5% when the surface is 
polished with 600 grit polishing paper. Even higher strengths 
could be expected with a finer polishing medium. Even if one 
were to accept a value in the range of 11,000 psi one must 
consider the effects of delayed failure under constant applied 
loads. Uner these loads, even if small, failure can occur after 
a given period of time has passed. This has been attributed to 
the mechanism of stress corrosion whereby water molecules from a 
humid environment enhance crack growth. Cracks continue to grow 
over a period of time until catastrophic failure occurs. 
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Therefore one needs a reliable method of predicting life 
times under applied loads. The fatigue parameters for making 
lifetime predictions must be obtained under the service 
environment, if at all possible and can be obtained from one of 
three types of experiments: crack velocity, static fatigue and 
dynamic fatigue. 

Crack velocity can be measured directly as a function of the 
stress intensity factor on specimens specifically designed for 
fracture mechanics experiments. These speicmens contain a 
macroscopic crack which allows accurate measurement of crack 
velocity and stress intensity factor. For purposes of failure 
prediciton crack velocity data is not as reliable as that from 
static or dynamic fatigue strength experiments because data from 
large, preformed cracks is not necessarily relevent to the 
propagation of microscopic cracks present in brittle materials. ^ 

Static fatiue data is generally obtained by measuring the 
time to failure of a large number of samples at several constant 
applied stresses. 

Dynamic fatigue is generated by measuring fracture strength 
as a function of stress rate. For this study dynamic fatigue 
was chosen to develop the time-to-failure equation necessary to a 
design diagram. This concept will be developed fully in the 
Technical Discussion section. 

Due to problems in obtaining samples in time to perform the 
actual experiments no data will be presented in the report. 
Instead attention will be paid to theoretical development of the 
design diagram and the proposed experimental procedure to obtain 
the necessary data for the diagram. 

Technical Discussion 


Initial Strength 

Modulus of Rupture (MOR) will be determined using the double 
ring bending strength test method. Figure 1 is a graphical 
representation of this technique. In this case a disc shaped 
specimen is loaded between two concentric rings. For limited 
forces, a tensile stress field will be set up in the central 
region of the convexly bent speicmen surface. Outside the load 
ring the radial and tangential stresses in the specimen decrease 
toward the edge, so that the possibility of fracture there is 
small. By increasing the force the tensile stress in the 
specimen center is increased at a constant rate until fracture 
occurs, with the expectation that the fracture is initiated in 
the region of the surface subjected to the maximum stress, 
underneath the load ring. 
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The major advantage of this technique, as compared to techniques 
such as 3 and 4 point bending, is the elimination of edge 
effects. 

To calculate the rate of increase of the bending stress in 
circular specimens the following is applicable: 

= /.OX ■ AF (i) 

A* S i 

Where : 

A f = Increase of test force, measured in the time interval 
near the fracture initiation 

AS- Specimen thickness 

1.08 = Numerical constant related to the load ring and 

specimen support ring diameters and Poisson's Ratio 

The bending strength ) belonging to the maximum force 

(Fmax) from equation (1) is: 

Kffrg = 1.08 Fmax (2) 


Where: 

- bending strength 
Fmax = greatest force 
S = specimen thickness 


Weibull Statistics 

In order to determine the distributions of strength Weibull - 
plots will be utilized. Fracture of brittle materials can be 
regarded as a statistical process, since the material has a 
distribution of flaw sizes. Weibull 5 described the mathematics 
for this process in terms of probability functions. His 
equations can take forms: 
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(3) 


F = 1 - exp (--3T )‘ nrrv 


or 


F = 1 - exp 

To 


(4) 


where : 


F = 




/m = 


probability that the material will fail at stress 
constants which describe the distribution 
constant (Weibull modulus) 


Equation (3) is the two-parameter model and represents an 
unbounded distribution. Equation (4) is the three-parameter 
model and represents a lower bounded distribution. 

A number of dynamic fatigue studies of glass and glass- 
ceramics have shown the two-parameters model to be acceptable in 
describing the flaw population and it's relation to material 
failure. • • Therefore, in this study equation (3) will be 
utilized. 

One can rearrange equation (3) to obtain a more workable 
equation. Starting with equation (3) 



F = 1 - exp ~j£ f* 

(3) 

- 1 + 

F = - exp (- T v*- 


1 - 

F = exp ( -ST 


(1 ” F)- 

exp ( Sl )'*" 






or 


1 = exp 

( JC )*“ 

(5) 

1 - F 

T t 

Taking In 

of both sides of 

(5) 
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and taking In again; 


ln ^ A/_Z__) =/yh/A}xr~/friJ^*TZ 

Where : 

"ST “ Failure stress (MOR) 

And: 

Z 17 - /?2 - O . jL, 

AJ 

Where : 

= Rank of sample in terms of strength 
Al - Number of samples tested 

Thus a plot of In ( In 1 ) versus ln^T" will yield a slope 

1-F 

of /?ft- . The slope indicates the distribution of strengths . As 

he slope (/flt) increases the distribution narrows, i.e. there is 
less scatter in the results. 

Lifetime Prediciton Equation 

Weibull statistics can be combined with data on static or 
dynamic fatigue to yield predictions of the lifetime of brittle 
materials under various environmental conditions. The 
derivation of this equation combines elements of fracture 
mechanics theory and Weibull statistics. 

The rate of crack cjrowth in a material can be expressed as a 
power law funciton: 9 ^ 


V 


Where : 


= AKj 


( 8 ) 


a = crack length 
Kj = stress intensity factor 
A, N = constant for a given material and environmental 
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The stress intensity factor can be expressed as 


Kj = TYE 

Where: 


(9) 


a = crack length 
CT = applied stress 

Y = crack geometry parameter which is constant for a 
given material 

Combining equations (8) and (9) to eliminate Kj yields: 




A y ^ 


( 10 ) 


For the case of constant applied stress ( ) : 


Ja = A 

d* 

or 

- a (id 


To obtain the material lifetime, this expression is 
integrated over the initial-to-f inal crack length and from the 
time the stress was applied = O) to the time of fracture 

( X = fp • 




(12) 
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6L 




+ i 


— A , /^2 -t I 




A) 'cl 


or 

/? ^y^T/ 


which yields: 


L _ 





a 


-C/O'Xjz)^ 



(13) 


Rearranging equation (9) : 



and substituting into equation (13): 



The stress intensity factor when fracture occurs, K lc , and 
the fracture stress is the stress at which fracture would 
occur in an inert enviornment unaffected by subcritical crack 
growth ( <ZJ2 ) 
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Substituting into equation (14) : 




Os) 


K 


-'Go-.i.) 

Io 


Generally ^ 


Therefore one can write: 

_= 2 . /05 y )"'* 

Solving for : 


(/£) 



-.A. a — s^xr-* (n) 

^ a M'Ia A 


Where : 


B = 


AY 2 


Taking In of of both sides of equation 


(17) 


!/ M 


^4' -t-(jO’'X) IsAQ'si — /[) //v \7X 


(18) 


Rearrangement of Equation (7) yields: 

IajQaj jir F J -t /7T^ l/o <^- a - /yy- /aj t - 


or 

IaJT = ^7 7^=)-+- N*- ^ 
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( 19 ) 


In an enviornment where no subcritical flaw grow occurs 
prior to fracture one can write equation (19) in terms of inert 
strength ( 


/a> fhQu-fe) -/ /s' VxJ 


Substitution of equation (20) into equation (18) for 1/J 
yields: 


!s> ^ (2D 


Where : 

represents the time required for a flaw to grow from an 
initial subcritical size to dimensions critical for catastrophic 
crack propagation leading to material failure. 

And: 

B and N are constants which characterize the subcritical flaw 
growth . 


Determination of B and N 

Evans 10 derived an expression relating the strength of a 
material ( ^7" ) to the stressing rate ( ) at constant 

temperature for the case of subcritical flaw growth as follows: 
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di 


or, 


dxr = <rd£ 

Dividing by da- ? 

r 

< vT 





z^r 


XT 


V 


Now, we know that, 
V « AK X 


Therefore we can write: 

dsn - 'j-da- 

. A/ 


AK-I 


(26) 


We also know: 

Kj = Y V45I 
or 

Kj = Y aJ '? a/ CL' 4// * 

Substituting this expression into equaiton (26) gives: 


a^T^ /j y v 


or 




( 22 ) 


(23) 


(24) 


(25) 
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( 27 ) 


Integrating over the extent of crack growth from an 
of zero to a final failure stress or obtains: 



% 

xr^dv = d- 

a y* 



a m da. 


which yields: 


AJ+i 


/? y+fr-A) j] 






O-'A* )/+ 


Since, 


a 


c. Ud 


the former can be dropped leading to: 


Md ( /) ysJ fa ) 

The intial flaw size can be taken as: 


af = k Ic /y n c 

or 

fa = j^L 

Substituting equation (31) into (30) for CZ* yields 


<% u+l 




intial stress 


(28) 


(29) 


(30) 


(31) 


( 32 ) 
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Can rewrite equation (32) as: 


y / m 

= B (N + l)<JZc V 


where : 


B = 


JL 


AYHV'X) Xi£~ x 

Taking In of both sides of equation (33) : 

(N + 1) 1a/ = 1w/b (N + 1) + (N - 2)^-JJ + 1 n^T 

or 

^ = In j^B (N + 1) + (N - 2 )\uj] + lA I O' 

Th - 7 


Rewriting Equation (34) : 


In ^ - 


1 N£b (N + lj + (N - 2) 1 a/<£ c 

JTTi ’~~ 


c) + XN'Y 

J N + 1 


or 


In ^ = Oo +<2|/w'T 


where: 

Oo = 

a, = 


IN [jB (N + 1)1 + (N - 2) 1 Njcc. 


Al+I 


N + 1 


(33) 

(34) 

(35) 

(36) 
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Thus, to obtain values for B and N, one fits fracture stress 
and stress rate data for samples tested in humid enviornment to a 
straight line g iven by equation (36), that is a plot 1 n^ versus 
1 amt .To ensure that the values of B and N are reliable one 
needs to test over a wide range of stress rates (i.e. over 
several orders of magnitude). The inert strength ( t TjL ) should 
be obtained at high enough stress rates to insure that crack 
growth is independent of environmental (i.e. no stress corrosion 
at the crack tip is occurring) . 

Design Diagrams 


Once the values of the constants B and N are determined, the 
time-to-failure for various values of failure probability (F) and 
applied stress (<Ta.) can be computed. Figure 2 shows an exmaple 
of a design diagram for polished ULE glass at F = 0.001, 0.01 and 
0.20. Comparing this to a diagram for unpolished ULE glass 
(Figure 3) one can note decrease in strength at a given time 
period (Ex. 10 yrs.) as compared to polished glass. This infers 
that at a given probability of failure the polished glass should 
withstand an applied stress for a longer period of time. For 
examples at F * 0.01 and an applied stress of 24 x 10° MPa the 
polished material will fail in one hour whereas the unpolished 
fails in one second. 

Experimental Approach 

The investigation will be carried out in three phases. 

Phase I- Residual Stress Investigation 

In order to determine the effects of residual stresses on 
the modulus of rupture of Zerodur, 10 samples with ground surfae 
finish will be annealed and compared with 10 duplicate unannealed 
samples using the double ring bend test method. All tests will 
be performed at a stress rate of 290 psi/sec, which is in 
accordance with test standard DIN 52292. 

Phase II - Generation of Design Diagram 

Weibull plots [1/0 ( ) vs. 1^^] will be made in order 

to determine the strength distributions for three types of sample 
preparation (AS cut, ground, and ground with acid etch). The 
etching technique is to be supplied by Schott Glass Co. MOR data 
will be obtained at four different stress rates, 3.5 x 10“ 3 
MPa/s, 2.7 x 10 “ 2 MPa/s, 2 MPa/s and 200 MPa/s. Thirty five 
samples will be tested. The final stress at each stress rate 
will be used to determine the inert strength. In this instance 
samples will be heated at 150°C for 24 hours in a dry N 2 
environment and put in a vacuum dessicator before testing. 
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Each sample will be exposed to flowing dry initrogen for fifteen 
minutes before testing, in a plexiglass chamber fitted to the 
testing machine. To insure that this stress rate leads to crack 
growth independent of environment, ten samples will be placed in 
distilled water for 24 hours and then tested at 200 MPa/s. The 
sample to be tested at the three lower stress rates will be 
stored at ambient conditions (T = 25°C, R.H. = 50%) for at least 
one week before testing. 

Using the time-to-failure equation developed in this paper, 
a time-to-failure diagram will be developed for each sample 
preparation. That is a separate diagram for the as cut, ground, 
and ground and acid etched samples. 

Phase III- Verficiation of Diagrams 

In order to verify the three time-to-failure diagrams, 
applied stress values which intersect the F = 0.005 probability 
curve corresponding to time-to-failures of 1 minute, 1 hour, and 
1 day will be selected and used to test 10 specimens each using 
the double ring test method. 
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